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Considering Daya Bay experimental value of sin ^13 we find tliat for standard interaction of neu- 
trinos witli matter the better discovery reach of CP violation is possible in short baseline neutrino 
oscillation experiment. However, this changes in presence of non-standard interactions (NSIs). For 
small non-standard interactions (small £0/3) long baseline is found to be more suitable, in general, 
for better discovery reach of CP violation. Even for large value of NSI long baseline is found to 
be better. We have also discussed discovery reach of hierarchy and the discovery reach of different 
NSIs for longer baseline. 



* zini@ctp-jamia.res.in 
t arnab@ctp-jamia.res.in 
rathin@ctp-jamia.res.in 



2 



I. INTRODUCTION 

Among the various neutrino oscihation parameters the values of two mixing angles 612 and 623 have been 
provided by experiments with certain accuracy. The magnitude of the mass square differences |Am|]^| and 
Am2i are also known but the sign of Arnl^ is still not known i.e. we still do not know exactly if neutrino 
masses follow normal heirarchy (NH) or inverted hierarchy (IH) . Apart from that we still do not know about 
the CP violating phase S. Earlier, the third mixing angle sin^ 2(?i3 was also not known accurately except its 
upper bound but recently the Daya Bay experiment [1] have predicted the value of this parameter with 5a 
confidence level. Now the only unknowns in the neutrino mixing matrix - so called PMNS matrix are CP 
violating phase 6 and the hierarchy of neutrino masses. 

In this work we consider neutrino superbeam (which mainly contains and P/j, ) coming from CERN 

travelling a baseline length of 2300 Km reaches Pyhasalmi (Finland) where a Liquid Argon detector is placed. 
We also consider another baseline of 130 Km with a superbeam source at CERN and Water Cherenkov 
detector placed at Frejus (France). We do a comparative study in the discovery potentials of the CP 
violating phase 6 for these two baselines in presence of both standard and non-standard interactions. 

This paper is organised as follows: In section II we discuss how the probability of oscillation for i/^ for 
short and long baseline differ when NSIs are taken into account. In section III we discuss the experimental 

setup and assumptions in doing the numerical analysis using GLoBES [2]. In section IV we discuss the 
discovery reaches of CP violation, hierarchy and NSIs. Finally in section V we conclude with remarks on 
the prospects of short and long baselines neutrino oscillation experiments. 



II. NEUTRINO OCILLATION PROBABILITIES WITH NSI 

In addition to the Standard Model (SM) Lagrangian density we consider the following non-standard 
fermion-neutrino interaction in matter defined by the Lagrangian: 

jCf^si = -2V2GFsi^^[h^Pf][uf,rLiya] (1) 

where P <E {L,R), L ~ ^^-^-^j R = ^^^2 ' ^ e,u,d and ej^^.^ is the deviation from SM interactions and 
can be called as the non-standard interactions (NSIs). A bound can be set to these NSI parameters [3, 4] 
which are dependent on specific models [5, 6] or are also model independent [7]. These NSI parameters can 
be reduced to the effective parameters and can be written as: 

^«/^ = E4f^ (2) 

where n/ is the fermion number density and Ug is the electron number density. These NSIs play significant 
role in the context of neutrino oscillation experiments and modify the interaction of neutrinos with matter 
and thus change the oscillation probability of different flavor of neutrinos. The NSIs could be present at the 
source of neutrinos, during the propagation of neutrinos and also during detection of neutrinos [8]. The NSI 
effects are expected to be smaller at the source and detector due to their stringent constraints [7]. We shall 
consider the NSI effect during the propagation of neutrinos only and in section IV in presenting results of 
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our numerical analysis wc shall consider the model independent allowed range of real values of different NSIs 
as mentioned in reference [7] for earth like matter. 



In vacuum, flavor eigenstates i/a may be related to mass eigenstates of neutrinos as 

\va>='^Uai\i'i> ; i = l,2,3, 



(3) 



where U is PMNS matrix [9] which depends on three mixing angles 612, 023 and dis and one CP violating 
phase 5. Although two more Majorana phases could be present in U but are not relevant for neurino 
oscillation experiments. The Hamiltonian due to standard (Hsm) and non-standard interactions (Hnsi) of 
neutrinos interacting with matter during propagation can be written in the flavor basis as: 



H = Hsm + Hnsi 
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In equations (5) and (6) 
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where m, is the mass of the i-th neutrino and A is considered due to the interaction of neutrinos with matter 
in SM, is the Fermi constant and Ug is the electron number density of matter . Sge, Sen > ^er; ^jur 
and Err arc considered due to the non-standard interaction (NSIs) of neutrinos with matter. In equation 
(6), ( * ) denotes complex conjugation. In general, the NSIs - e^^, and £ij_t could be complex. However, 
in our numerical analysis wc have considered those to be real. 

For understanding qualitatively which baseline is suitable for discovery reach of CP violation and hierarchy 
in presence of standard and non-standrad interactions we present below the oscillation probability P^^^^^ 
for both short and long baseline, — )• Vg oscillation channel is particularly sensitive to CP violation. In our 
numerical analysis however, we have considered other channels also like j/^ — > t'p. To get these expressions 
of probability we have followed the perturbation method adopted in references [8, 10-12]. At first, we shall 
present the oscillation probability upto order as the CP violating phase S appears at this order for only 
SM interactions as well as for small NSI of the order of a. For short baseline of 130 Km keeping the matter 
effect parameter A at the 0{a) and considering NSI parameters Ec/S of the order of a one obtains 
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For 2300 Km the matter effect is large. Keeping A in the leading order of perturbation and considering NSI 
parameters of the order of a one obtains 
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For CP violation there is difference of probability in the neutrino oscillation and probability of antineutrino 
oscillation. One can relate the oscillation probabilities for antineutrinos to those probabilities given for 
neutrinos above by the following relation: 



PocS = Pa^iScp -Sep, A -)■ -A). 



(11) 



In addition, we also have to replace with their complex conjugates, in order to deduce the oscillation 
probability for the antineutrino, if one considers non-standard interaction during propagation. 

To estimate the order of magnitude of different terms in the above two oscillation probabilities we shall 
consider A ^ a for short baseline and A ^ 0.5 for longer baseline and following Daya Bay result we 
shall consider sin(?i3 ^ ^/a. For only SM interactions, (i.e Sap 0) in above expressions of oscillation 
probabilities one finds that for both short and long baseline the 6 dependence occurs at order of a"^/^. 
This order of dependence with 6 remains same in the difference of neutrino oscillation probabilities and 
antineutrino oscillation probabilities (represented by AP later). However, the S independent part in AP 
(which could mimic CP violation) for short baseline is at order but for long baseline this is at order a. 
This happens due to matter effect through A for SM as can be seen from above expressions. For this reason 
the discovery reach of CP violation is better in short baseline than that in long baseline. However, when 
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NSIs arc also taken into account one can sec that for longer baseline further S dependence in AP could occur 
at the order of a^^^ through 02 and 03 containing terms in (9) for NSIs of the order of a. We have checked 
that for slightly higher NSIs of the order of ^/a using perturbation method the same 5 dependent terms in 
AP appears with 02 and 03 in the oscillation probability for long baseline as given in (9) and this slightly 
higher NSI makes these terms at the order of a which could compete with the S independent part (which 
could mimic CP violation) in AP for long baseline as that is also at the order of a. This improvement of 5 
dependent part over independent part for long baseline does not happen for short baseline. First of all for 
NSI of the order of a oscillation probability in (8) for short baseline is independent of NSIs. Even for slightly 
higher NSIs of order we have checked although NSIs like 02 and as enter into the only S dependent 
part of AP for short baseline but that is at the order of a^. So not much improvement in the discovery 
reach of CP violation is expected for short baseline for NSIs which are present in 02 and 03. So presence of 
slightly higher NSIs of order ^/a present in 02 and as improves the discovery reach of CP violation for longer 
baseline in comparison to the short baseline. As 02 and 03 contains NSIs like and Egr it is expected that 
in presence of these NSIs the long baseline could provide a better discovery reach for CP violation. 

For hierarchy discovery the difference in the oscillation probabiliticis dnc to two different hierachies mat- 
ters. If we change hierarchy then in above oscillation probabilities the following transformations are to be 
considered: 

Amli -Amgi ; A -A ; a -a 

For SM only the above differences of oscillation probabilities between two different hierarchies (APfj) for 
longer baseline is of the order a due to matter interaction through A. However, for shorter baseline the AP^ 
is of the order a^/^. So longer baseline is more suitable for discovery of hierarchy. This remains unaltered 
even after including NSIs. 



III. ANALYSIS AND EXPERIMENTAL SETUP 

In this work for the numerical simulation we consider two set-ups: (a) A Superbeam setup which originates 
in CERN and reaches a 100 kt Liquid Argon detector placed at a distance of 2300 Km at Pyhasalmi (Finland) 
(b) A Superbeam setup originating in CERN and reaching a 500 Kt Water Cherenkov detector placed at a 
distance of 130 Km at Frejus (Prance). 

We consider a flux with mean neutrino energy around 5 GeV and 3 x 10^^ protons on target per year. We 
consider the same flux as in [13]. For set-up(a) we consider a power of 0.8 MW per year. In doing the analysis 

we consider a signal efficiency of 90% in both the appearance and disappearance channels. As backgroiuids 
we consider a 0.5% neutral current events, 1% of misidentified and the complete intrinsic contamination 
of the beam. The background rejection efficiencies were assumed to be constant over the energy window of 
0.5 to 10 GeV. The calibration error has been considered to be 2% for disappearance channel and 0.01% 
for Ue appearance channel. The Gaussian energy resolution is considered to be 150 MeV for electrons and 
positrons and for muons. The correlation between the visible energy of background NC events and 

the neutrino energy is implemented by migration matrices which has been provided by L. Whitehead [14]. 

In the case of set-up(b) we have considered a beam intensity of 4 MW and a systemetics on signal and 
background as 2% based on ref[15] . A neutrino energy window of 0.5 to 10 GeV has been considered for both 
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the set-ups and both v and v beams have been used simultaneously for a time period of 5 yrs for neutrino 
and 5 yrs for antineutrino in case of set up (a) 2300 km and 2 yrs for neutrino and 8 yrs for antineutrino in 
case of set up (b) 130 km. 

Following ref [13] we consider the true values of the neutrino oscillation parameters as | Am|]^ | = 2.45 x 10"*^ 
eV^, Amii = 7.64 X 10"^ eV^, 6*13 = 9°, 6*12 = 34.2° and 6*23 = 45°. Also in calculating the priors we consider 
an error of 3% on ^12, 0.005 on sin^20i3, 8% on 623, 4% on [Amlil and 2.5% on Am2i. Also we consider 
an error of 2% on matter density and a systematic uncertainties of 5 and 10% on signal and background 
channels. In doing the whole analysis we have used GLoBES software [2]. In this work we have done a 
comparative study of the two set-ups (a) and (b) in finding the discovery reach of the CP violating phase 5 
depending on different NSI parameters and the hierarchies. For the longer baseline we have also discussed 
about discovery reach of hierarchy and NSI. 



IV. RESULTS 
A. Discovery of CP violation 
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FIG. 1: CP violating phase 5 discovery for two different baselines 130 Km and 2300 Km. The right hand panel is 
for IH and the left hand panel is for NH . 



In figure 1 we have shown the discovery reach oi CP violation for SM interactions of neutrinos with matter. 
For 130 Km baseline the CP violation could be discovered over about 67% of the possible 6 values for normal 
hierarchy and this is about 68% for inverted hierarchy whereas for 2300 Km baseline these values are about 
46% and 47% respectively. The discovery reach for longer baseline of 2300 Km was shown earlier by Coloma 
et al [13]. So with only SM the short baseline like 130 Km seems to be better for good discovery reach of CP 
violation. However, this could change depending on the type of NSIs taken into account as discussed below. 
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FIG. 2: CP violating phase 6 discovery for two different baselines 130 Km and 2300 Km at 5cr considering NSIs Eee, 
een and Ser- The right hand panel is for IH and the left hand panel is for NH . 



In figure 2 we have compared the discovery reach of CP violation for 130 Km baseline and 2300 Km 
baseline in presence of NSIs like Sge, Ce/n and Eer • For Sge in the range of about - 2.5 to 0.5 the 2300 Km 
baseline is found to be better for the discovery reach which could be possible for about 49 % of the possible 
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FIG. 3: CP violating phase 5 discovery for two different baselines 130 Km and 2300 Km at 5cr considering NSIs e^^^, 
e^T and err- The right hand panel is for IH and the left hand panel is for NH . 



i5 values whereas 130 Km is not so good in this range of NSI See- For somewhat larger NSI aroimd 4 or -4 
the short baseline of 130 Km is found to be better for discovery reach which could be possible for about 
36% of the possible 6 values. It is interesting to note that for NSI - Sen - the discovery reach is consistently 
better for 2300 Km than that for 130 Km and could be discovered over 90% of the allowed 5 values. For Ser 
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in 2300 Km CP violation could be discovered over about 70% of the 5 values over the entire allowed region 
of £ex- 130 Km baseline is found to be better for very high s^t values around 3 and -3 only. 

In figure 3 we have compared the discovery reach of CP violation for 130 Km baseline and 2300 Km 
baseline in presence of NSIs like e^^jE/ur and Ser ■ For particularly and e^r, 130 Km baseline is found 
to be better for CP discovery which could be possible over 70 % of the 6 values for both the hierarchies 
throughout the entire allowed region of these NSIs. This is about 43 % for 2300 Km baseline. For Ett in 
2300 Km baseline the discovery reach of CP violation could be possible over 50 % of (5 values for small NSIs. 
For 130 Km for large NSIs only good discovery reach of about 61% could be possible. So it is found that 
except e^^ and e^^ for other NSIs for their small values the discovery reach for CP violation is better for 
longer baseline than the shorter baseline. 



B. Discovery of hierarchy 
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FIG. 4: Discovery reach of hierarchy for two different baselines 130 Km and 2300 Km for SM interactions only. The 
right hand panel is for IH and the left hand panel is for NH . 



In figure 4 it is seen that with only SM interactions the discovery reach of hierarchy could be possible at 
above 5a confidence level for 2300 Km baseline. It was shown earlier by Coloma et al [13]. The 130 Km 
baseline is found to be not so suitable for hierarchy as the discovery could be possible at only below 3(7. 
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FIG. 5: Discovery reach of hierarchy with NSIs for a baseline of 2300 Km at 3a. 



In figure 5 the discovery reach of hierarchy is shown in the d NSI plane. Although for entire negative 
allowed values of See hierachy could be discovered at 3 cr confidence level but for positive values it could be 
possible for lesser than 2. The discovery reach is possible for the entire allowed region of NSIs particularly 
for Egfj,, Ser and £^^. For other NSIs like e^x and Srr discovery is possible for almost entire allowed region 
except few S and NSI values found in small enclosed areas. For Srr near to about 1 the discovery is difficult. 
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C. Discovery of NSI 




En, (true) 

FIG. 6: Discovery reach of NSIs for both the hierarchies for a baseline of 2300 Km at 5a. 



In figure(6) we have shown the discovery reach of difFcrent NSIs with 6 for both the hierarchies at 5a 
confidence leveL The NSI discovery reach in the figures are the area bounded by the sofid (IH) and dotted 
(NH) lines for set-up(a) i.e. basehne is 2300 km. However, for set-up(b) i.e. baseUne 130 Km (not shown 
in figure) we are not able to obtain any statistically significant bound on the NSIs in the total allowed 
range of NSIs. From the figure the range (after approximately averaging for different delta values and 
also for different hierarchies) over which the different NSIs could be discovered for set-up (a) are given as 
-0.32 < See < 0.32, -0.038 < < 0.036, -0.07 < Ser < 0.075, -0.26 < e^^^^ < 0.5, -0.016 < < 0.016 
and -0.3 < e^r < 0.26. 
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V. CONCLUSION 



In this work we have considered two experimental set-ups - one with a basehne of 2300 Km directed 
towards a 100 Kt Liquid argon detector and the other with a basehne of 130 Km directed towards a 500 Kt 
Water Cherenkov detector. Considering the recently obtained value of sin ^13 from Daya Bay experiment 
it is found that in presence of only SM interactions the discovery reach of CP violation is better for short 
baseline like 130 Km than long baseline like 2300 Km. In presence of NSIs e^^ and Sf^r for their entire 
allowed region the short baseline is found to be more suitable for discovery of CP violation. However, for 
£e/i over its entire allowed region and also for Ser over significant part of the allowed region in the 6 NSI 
plane the longer baseline is found to be better. For other NSIs in general for smaller values the longer 
baseline is better whereas for larger NSI values the shorter baseline is found to be better. In presence of 
only SM interactions the discovery reach of hierarchy is far superior for longer baseline than that for short 
baseline. In presence of NSI in long baseline we find that the discovery reach of hierarchy is possible for the 
entire allowed region of NSIs particularly for Eeju, ^er and 6^^. However, for allowed region of £ee < 2 only 
it is possible to find the good discovery reach at 3(t. For other NSIs like and Err discovery is possible 
for almost entire allowed region except for few 6 and NSI values. The discovery reach of various NSIs are 
possible in the longer baseline for NSI values greater than the order of a. 

The good discovery reach of CP violation crucially depends on the relative contribution of the 6 dependent 
terms in the probability with respect to the d independent part as discussed in section II. In the longer baseline 
the NSIs Sen and Eer is coupled with S in the oscillation probability and as such increases the S dependent 
contribution and provides the possibility of good discovery reach of CP violation. This does not happen for 
very large NSI values and the shorter baseline is found to be better for CP violation discovery. One may 
note that the expression of oscillation probability which we have presented in section II are not valid for 
such high NSI values. Considering the possibility of the presence of NSIs in nature it seems both short and 
long baseline should be considered particularly for discovery of CP violation in the leptonic sector through 
neutrino oscillation experiment. 
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